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Abstract 

Background: Early-life exposure to environmental microbial agents may be asso- 
ciated with the development of allergies. The aim of the study was to identify bet- 
ter ways to characterize microbial exposure as a predictor of respiratory 
symptoms and allergies. 

Methods: A birth cohort of 410 children was followed up until 6 years of age. 
Bacterial endotoxin, 3-hydroxy fatty acids, A r -acetyl-muramic acid, fungal extra- 
cellular polysaccharides (EPS) from Penicillium and Aspergillus spp., (3-D-glucan, 
ergosterol, and bacterial or fungal quantitative polymerase chain reactions 
(qPCRs) were analyzed from dust samples collected at 2 months of age. Asthma, 
wheezing, cough, and atopic dermatitis were assessed using repeated question- 
naires. Specific IgEs were determined at the age of 1 and 6 years. 
Results: Only few associations were found between single microbial markers and 
the studied outcomes. In contrast, a score for the total quantity of microbial expo- 
sure, that is, sum of indicators for fungi (ergosterol), Gram-positive (muramic 
acid) bacteria, and Gram-negative (endotoxin) bacteria, was significantly 
(inverted-U shape) associated with asthma incidence (P < 0.001): the highest risk 
was found at medium levels (adjusted odds ratio (aOR) 2.24, 95% confidence 
interval (95% CI) 0.87-5.75 for 3rd quintile) and the lowest risk at the highest 
level (aOR 0.34, 95% CI 0.09-1.36 for 5th quintile). The microbial diversity score, 
that is, sum of detected qPCRs, was inversely associated with risk of wheezing 
and was significantly (inverted-U shape) associated with sensitization to inhalant 
allergens. 

Conclusion: Score for quantity of microbial exposure predicted asthma better 
than single microbial markers independently of microbial diversity and amount of 
dust. Better indicators of total quantity and diversity of microbial exposure are 
needed in studies on the development of asthma. 



Microbial exposure may contribute to the development of 
allergic diseases, but the role of different agents and mecha- 
nisms are still unknown (1). Most epidemiological studies 
with microbial exposure assessment among children have 
been restricted to studies on endotoxin (2) or have used only 



a limited set of microbial markers (3-7). These microbial 
markers reflect typically only a narrow part of the total 
microbial exposure. Recently, some studies have applied 
combined indicators of microbial exposure (8, 9), like diver- 
sity of microbes, which appears to have a protective effect on 
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the development of asthma (8). Thus, there is a clear need to 
use combined indicators and other new approaches to 
describe exposure to environmental microbes. 

The objective of this study was therefore to study prospec- 
tively the associations between several markers of environmen- 
tal microbial exposure — including their combinations — and 
the development of doctor-diagnosed asthma, respiratory 
symptoms, atopic dermatitis, and sensitization up to the age of 
6 years. 

Methods 

The study population consisted of children born in Eastern 
and Middle Finland: the first half of the study population 
(# = 214) belongs to a European birth cohort (PASTURE) 

(10) among farmers and nonfarmers, while the second half of 
the cohort is an extension to unselected children (N = 228) 

(11) . Pregnant women were recruited, whose children were 
born between September 2002 and May 2005. The study pro- 
tocol was approved by a local ethics committee in Finland 
(11). A written informed consent was obtained from the 
parents. 

Follow-up 

The first questionnaire was administered during the third 
trimester of pregnancy (11). The follow-up questionnaires 
(at the age of 2, 12, 18, and 24 months, and thereafter 
annually), except the questionnaire at 2 months, enquired 
about any wheezing, other respiratory symptoms apart from 
cold (wheezing, cough, nocturnal cough), and confounders 
for the time period after the preceding questionnaire. 
'Asthma ever' was defined as first parent-reported doctor- 
diagnosed asthma and/or second diagnoses of asthmatic (or 
obstructive) bronchitis. 'Current asthma' was defined as 
those asthma ever cases who also reported the use of 
asthma medication and/or wheezing in past 12 months in 
the 6-year follow-up. 'Atopic dermatitis ever' was defined as 
first parent-reported doctor-diagnosed atopic dermatitis dur- 
ing the follow-up. 

House dust samples 

The protocols for dust collection are explained in detail in 
the Supplementary material. Briefly, in the Finnish PAS- 
TURE study, fieldworkers took two samples from living 
room floors for the present report. The first samples were 
processed without sieving at Utrecht University, where they 
were extracted and analyzed for endotoxin (as a marker of 
Gram-negative bacteria) and EPS as described previously 
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(12) , and the stored residue after extraction was later ana- 
lyzed for P(l,3)-glucans (as a marker of fungal biomass) by 
heat extraction and inhibition enzyme immunoassay (EIA) 

(13) . The second samples were processed and analyzed at 
National Institute for Health and Welfare (THL), where they 
were sieved to remove larger particles, dried, split and then 
stored, and finally analyzed for muramic acid (as a marker of 
Gram-positive bacteria), ergosterol (as a marker of fungal 
biomass), and lipopolysaccharide (LPSi 0: o-i6:o) ( as a marker 
of Gram-negative bacteria) by gas chromatography tandem 
mass spectrometry (GC-MS-MS) (14) and by quantitative 
polymerase chain reaction (qPCR) for microbial DNA 
(15-17). 

In the extended cohort, results from only one dust sample 
from the living room floor collected by the families are used 
in the present report. The sample was sieved and dried at 
THL, then split in aliquots for various analyses: One aliquot 
was sent to Utrecht where it was extracted and analyzed for 
endotoxin, EPS, and (3(l,3)-glucans (13), and other aliquots 
were analyzed for the chemical markers (muramic acid, 
ergosterol, LPSi 0: o-i6:o) (14) and microbial DNA (15-17). 

Quantitative polymerase chain reaction (qPCR) (15-17) 
was used to determine microbes that have been suggested to 
be common in moisture-damaged buildings: two genera of 
Gram-positive bacteria, Mycobacterium spp. and Streptomy- 
ces spp., and six fungal species, genera, or groups: the assay 
group for Aspergillus fumigatus/Neosartorya fiseheri (AspNeo 
group); Cladosporium spp.; the combined assay group for 
Penicillium spp., Aspergillus spp., and Paecilomyces variotii 
(PenAsp group); Stachybotrys chartarum; Trichoderma viridej 
atroviride/koningii (Trichoderma viride group); and Wallemia 
sebi. 

Exposure to all markers was expressed as loads, that is, the 
amount per m sampling surface. As the laboratory proce- 
dures determine concentrations (agents per mg dust), loads 
are calculated by multiplying the concentrations with the mea- 
sured dust weight per m 2 . However, as endotoxin, EPS, and 
glucans in the PASTURE samples were determined in nonsi- 
eved dust, while all other analyses were performed with sieved 
dust, we calculated loads for the former parameters by multi- 
plying with the load of sieved dust in the parallel sample taken 
for the other analyses. Weights for these parallel sieved and 
nonsieved dust samples in the Finnish PASTURE study 
showed a fair correlation (Spearman r = 0.7), but with a sys- 
tematic difference: as expected, sieved dust samples contained 
on average 60% (95% CI 57, 64) less dust. Mean concentra- 
tions in sieved and nonsieved dust from nonfarm homes were, 
however, very similar, which suggests that endotoxin, EPS, or 
glucans are not preferentially bound to the smaller (sieved) or 
larger (removed by sieving) particle fractions. 



Abbreviations 

EPS, fungal extracellular polysaccharide antigens of Penicillium and Aspergillus species; GC-MS-MS, gas chromatography tandem mass 
spectrometry; LPS, lipopolysaccharide; LPSio : o-:i6:o:o. amount (mol) of LPS in each sample. The number of mol of 3-hydroxy fatty acids of 
C 10:0 to C 16:0 carbon chain length divided by 4; qPCR, quantitative polymerase chain reaction; PASTURE, Protection against Allergy STUdy in 
Rural Environments; OR, odds ratio; aOR, adjusted odds ratio; slgE, specific immunoglobulin E; kU/L, kilo units per liter; GEE, general 
estimating equations; EIA, inhibition enzyme immunoassay. 
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Blood sampling 

Venous blood samples collected at one (N = 374) and at 
6 years of age (N = 304) were analyzed for specific immuno- 
globulin E (slgE) to 13 inhalant (Dermatophagoides pteronys- 
sinus, D. farinae, alder, birch, European hazel, grass pollen 
mixture, rye, mugwort, plantain, cat. horse and dog dander, 
and Alternaria alternate!) and 6 food allergens (hen's egg, 
cow's milk, peanut, hazelnut, carrot, and wheat) (Mediwiss 
Analytic, Moers, Germany) (18). The cutoff levels for slgE 
concentrations were 0.35 kU/1 at one and 0.70 kU/1 at 
6 years of age. The higher cutoff was used at the age of 
6 years due to high prevalence of sensitization (Table 1). 

Statistical analyses 

Generalized estimating equations (GEE) with an exchange- 
able correlation structure to account for correlation between 
repeated measures within subjects were used to determine 
associations between microbial exposure and repeated mea- 
sures of parent-reported wheezing and cough at different 
ages. Survival analyses (discrete-time hazard models) were 
used in analyzing asthma ever, current asthma, and atopic 



dermatitis ever. Logistic regression was used to analyze sensi- 
tization to inhalant allergens. The results are presented as 
adjusted odds ratios (aORs) and their 95% confidence inter- 
vals (95% CI). 

In the analysis of single microbial markers, the amounts of 
dust and microbial markers were generally divided into three 
categories using tertiles as cutoffs. For Wallemia sebi, the 
lowest category consisted of levels below the limit of detec- 
tion (37.4%) and the rest of the values were evenly divided 
into two groups (medium and high). Staehybotrys chartarum 
and AspNeo group were dichotomized at the detection level 
due to high numbers of nondetected values. Microbial diver- 
sity score was defined as the sum of all the detected qPCR 
markers (range 0-8). Due to a low number of observations 
with low and high diversity, the variable was categorized into 
four classes: 0-4, 5, 6, and 7-8. 

Total microbial quantity score was calculated as the sum 
of three markers for the three different microbial groups, that 
is, for fungi, Gram-negative bacteria, and Gram-positive bac- 
teria (one marker per group). Single markers were first 
divided into five groups using quintiles as cutoffs with scores 
0, 1, 2, 3, and 4 from lowest to highest and then summed. 



Table 1 Incidence of doctor-diagnosed asthma and atopic dermatitis and the point prevalence of respiratory symptoms during the first 
6 years and prevalence of sensitization to inhalant allergens at the age of 1 and 6 years 





At 1 


year 




At 1 


.5 years 


At 2 years 




At 3 years 




At 4 


years 




At 5 


years 




At 6 


years 




N 


n 


% 


N 


n 


% 


N 


n 


% 


N 


n 


% 


N 


n 


% 


N 


n 


% 


N 


n 


% 


Asthma ever* 


389 


6 


1.5 


353 


9 


2.6 


332 


14 


4.2 


291 


15 


5.2 


271 


9 


3.3 


257 


6 


2.3 


247 


3 


1.2 


Current 


353 


3 


0.9 


331 


7 


2.1 


315 


8 


2.5 


290 


7 


2.4 


281 


3 


1.1 


274 


4 


1.5 


270 


3 


1.1 


asthmaf 












































Wheezing 












































Apart from 


391 


41 


10.5 


358 


24 


6.7 


358 


25 


7.0 


338 


23 


6.8 


358 


29 


8.1 


346 


22 


6.4 


357 


25 


7.0 


cold 












































Any 


391 


91 


23.3 


358 


58 


16.2 


358 


67 


18.7 


339 


65 


19.2 


358 


68 


19.0 


348 


46 


13.2 


357 


51 


14.3 


Cough 












































Apart from 


391 


98 


25.1 


358 


66 


18.4 


358 


85 


23.7 


338 


80 


23.7 


357 


91 


25.5 


348 


87 


25.0 


357 


94 


26.3 


cold 












































Nocturnal 


386 


18 


4.7 


358 


47 


13.1 


358 


76 


21.2 


338 


63 


18.6 


357 


78 


21.9 


348 


59 


17.0 


357 


79 


22.1 


Atopic 


390 


61 


15.6 


300 


24 


8.0 


265 


12 


4.5 


230 


14 


6.1 


21 I 


9 


4.3 


199 


7 


3.5 


188 


2 


1.1 


dermatitis evert 












































Sensitization 


374 


78 


20.9 
































304 


115 


37.8 


to inhalant 












































allergens§ 













































N = in respirtory symptoms: total number of children with information on given symptom at each follow-up; in sensitization: total number of 
children with IgE results at a given time point; in asthma ever: total number of children at each follow-up who has information on asthma 
ever at each follow-up and has not censored from the follow-up due to onset of asthma ever; in current asthma: total number of children at 
each follow-up who has information on asthma ever and asthma medication/any wheezing at each follow-up and has not censored from the 
follow-up due to onset of current asthma; in atopic dermatitis ever: total number of children at each follow-up who has information on atopic 
dermatitis at each follow-up and has not censored from the follow-up due to onset of atopic dermatitis, n = the number of children with the 
presence of wheezing/cough symptom or with positive IgE result, or the number of newly diagnosed asthma/current asthma/atopic dermati- 
tis cases at each follow-up. 

* Doctor-diagnosed asthma ever is defined as doctor-diagnosed asthma at least once and/or doctor-diagnosed asthmatic bronchitis more than 
once during the 6-year follow-up. 

f Current asthma is defined as asthma ever with asthma medication and/or wheezing symptom past 12 months at the age of 6 years. 
^Atopic dermatitis ever is defined as docotor-diagnosed astopic dermatitis during the follow-up. 

§Atopic sensitization is defined as any inhalant slgE >0.35 kU/l at the age of 1 year and >0.70 kU/l at the age of 6 years. 
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Environmental microbial exposure and asthma 



All possible combinations with measured microbial markers 
were used to create four different total quantity scores for 
microbial exposure: (1) ergosterol, endotoxin, and muramic 
acid; (2) fS-D-glucan, endotoxin, and muramic acid; (3) ergos- 
terol, LPS 10: o-i6:o, and muramic acid; and (4) (3-D-glucan, 
LPS 10: o-i6:o, and muramic acid. The created variables, with 
ranges from 0 to 12, were divided into five groups using 
quintiles as cutoffs for statistical analyses. Linear and qua- 
dratic trends in the associations between microbial diversity 
or total quantity scores and health outcomes were tested 
using polynomial contrasts (19). 

All models were adjusted for study cohort, farming, and 
well-known risk factors for asthma (maternal history of aller- 
gic diseases, gender, number of older siblings, smoking dur- 
ing pregnancy). The models of total microbial quantity score 
(ergosterol, endotoxin, and muramic acid) and asthma ever, 
atopic dermatitis ever, and sensitization to inhalant allergen 
were carefully tested for 20 additional confounding factors, 
which have been described earlier (12). If a confounder chan- 
ged the estimates of the total microbial quantity score by 
more than 10%, all the models of the outcome were adjusted 
for the factor. All the models of respiratory symptoms and 



current asthma were adjusted for the same confounding fac- 
tors as in the asthma ever models. Given an exposure preva- 
lence of 25%, power calculations ((5 = 80%, a = 0.05) show 
that with these data, we can detect an odds ratio of 2.52 for 
asthma ever with a population prevalence of 14%. The data 
were analyzed using SAS 9.2 for Windows (SAS Inc., Cary, 
NC, USA). 

Results 

The present analyses included all those 410 children with 
data on at least one microbial marker available in any of the 
quantity scores. Of those, 341 (83%) had the results of all 
microbial determinations. During the 6-year follow-up, 62 of 
389 children with data on asthma developed asthma, and 
of those children. 35 still had current asthma at the age of 
6 years (Table 1). Distributions of all microbial markers are 
shown in Table 2. 

Most correlations between the amount of house dust and 
single microbial markers as well as intercorrelations between 
two microbial markers together were between 0.30 and 0.85 
(Table SI). Total microbial quantity scores were correlated 



Table 2 Levels of microbial markers 



Percentiles 





N 


<LOD* 


(%) 


25th 


Median 


75th 


95th 


Amount of dust (mg/m 2 )t 


410 


0 




173 


311 


587 


1332 


Gram-positive bacteria or marker of Gram-positive bacteria 












Mycobacterium (cells/m 2 ) 


399 


21 


(5.3) 


186 857 


641 427 


1 753 136 


7 041 045 


Streptomyces (cells/m 2 ) 


399 


7 


(1.8) 


3923 


13 089 


45 939 


462 289 


Muramic acid (ng/m 2 ) 


386 


7 


(1.8) 


3122 


6616 


13 713 


34 852 


Marker of Gram-negative bacteria 
















LPS 10:0 i6:o (nmol/m 2 ) 


398 


2 


(0.5) 


3 


8 


15 


46 


Endotoxin (EU/m 2 ) 


372 


0 




2668 


6253 


13 485 


66 569 


Fungi or marker of fungi 
















PenAsp group (cells/m 2 ) 


399 


3 


(0.8) 


13 580 


59 405 


233 616 


1 282 236 


Trichoderma viride group (cells/m 2 ) 


400 


62 


(15.5) 


340 


2028 


9906 


153 816 


Wallemia sebi (cells/m 2 ) 


400 


152 


(38.0) 


0 


883 


11 541 


1 52 041 


Cladosporiwn (cells/m 2 ) 


399 


50 


(12.5) 


3652 


25 683 


1 04 690 


982 138 


AspNeo group (cells/m 2 ) 


400 


320 


(80.0) 


0 


0 


0 


14 254 


Stachybotrys chartarum (cells/m 2 ) 


400 


390 


(97.5) 


0 


0 


0 


0 


EPS (EPSU/m 2 ) 


374 


15 


(2.1) 


7079 


15 312 


38 025 


126 107 


p-D-glucan (ng/m 2 ) 


376 


0 




511 


970 


1818 


4066 


Ergosterol (ng/m 2 ) 


401 


0 




1473 


3015 


7040 


21 978 


Quantity scores 
















Ergosterol, endotoxin, 


347 






3 


6 


9 


12 


and muramic acid 
















p-D-glucan, endotoxin, 


352 






3 


6 


9 


12 


and muramic acid 
















Ergosterol, LPS 10: o-i6:o. 


380 






3 


6 


9 


12 


and muramic acid 
















p-D-glucan, LPS 10: o 16:0. 


354 






3 


7 


9 


12 


and muramic acid 
















Diversity scorej 


399 






5 


6 


6 


7 



*LOD = the level of detection, percentage of samples under detection level from analyzed samples (%). 

fThe amount of dust refers to dust samples used for analyses of qPCRs, ergosterol, muramic acid, and LPS (sieved). 

JJThe number of detected 6 fungal or 2 bacterial qPCRs. 
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Amount of dust 

2nd fertile vs. 1st fertile 
3rd tertile vs. 1st tertile 

Mycobacterium (qPCR) 
2nd tertilfl vs. 1st tertile 
3rd tertile vs. 1st tertile 

Sireptomyces (qPCR) 
2nd tertile vs. 1st tertile 
3rd tertile VS. 1st tertile 

Muramic acid 
2nd tertile vs. 1st tertile 
3rd tertile vs. 1st tertile 

LPS10:0-16:0 
2nd tertile vs. 1st tertile 
3rd tertile vs. 1st tertile 

Endotoxin 
2nd tertile vs. 1st tertile 
3rd tertile vs. 1st tertile 



Asthma ever 


Any wheezing 


Nocturnal cough 


Sensitization tn inhalant allergen at ii years 


















1 •- 


1 


1 • 




1 • 


1 1 


» 1 


1 • 


^ 


1 • 


1 — • 


1 


1 • 


1 


1 • — 


— i 


1 •— 


— 1 


1 • 


1 


1 •■ 


1 


















1 < 


i 1 


1 •- 


— 1 


t- 


— • 1 


1 • 


H 


1 




1 • 1 


1 < 


l 1 


1 • 1 


1 

| * 


• 1 

H 


1 * 

1 • 


1 

— 1 


1 •— 

1 • 


H 

1 


1 

1 


« 1 

« 1 


















1 — 


— • 1 


1 • 


— 1 


1 • 


i 1 


1 • 


H 


1 • 


— 1 


1 • 




1 • 1 


1 • 




h 


♦ 1 


1 


-» 1 


>- 


— < 


1 


— • 1 


1 • 


1 


1 1 


» 1 


1 


• 1 


1 ' 


» 1 


l i 

0.25 0.5 


ii ii 

2 4 0.25 0,5 


ii il 

2 4 0.25 0.5 


II II 
2 4 0.25 0,5 


1 i 

2 4 



Odds ratios with 95% CI 



PenAsp grouplqlTR) 
2nd tertile vs. Isl tertile 
3rd tertile vs. 1st tertile 

Triehodemtaviride group (qPCR) 
2nd tertile vs. 1st tertile 
3rd tertile vs 1st tertile 

Wallemui sebi (qPCH) 
2nd tertile vs. 1st lertilc 
3rd tertile vs. 1st lertilc 

Cladosporium (qPCR) 
2nd tertile vs. 1st tertile 
3rd tertile vs. 1st tertile 

AspNeo group (qPCR) 
detected (20%) vs. not detected (SOW) 

Slachybotrys chartorum (qPCR) 
detected {25%) vs. not detected (97.3%) 

EPS 

2nd tertile vs. 1st tertile 
3rd tertile vs. 1st tertile 

G-D-gJucan 
2nd tertile vs. 1st tertile 
3rd tertile vs. 1st tertile 

Ergoslerol 
2nd tertile vs. 1st tertile 
3rd tertile vs. 1st tertile 



I • 1 

I • 1 



I • 1 



I • 1 

I • 1 



I • 1 



Any wheezine. 



I • 1 

I • 1 



I • 1 



I • 1 

I • 1 



2 4 8 0.25 0.5 1 
Odds ratios with 95% CI 



0.25 0.5 1 2 4 



Figure 1 Adjusted associations between amount of dust and bac- 
terial (A) and fungal (B) exposure at the age of 2 months and 
asthma and respiratory symptoms during the first 6 years of life 
and sensitization to inhalant allergen at the age of 6 years. Models 
are adjusted for study cohort, farming status, gender, maternal his- 



tory of allergic diseases (hay fever, atopic dermatitis, and/or 
asthma), smoking during pregnancy, the number of older siblings, 
and living area of the home. Models of sensitization to inhalant 
allergens (cut point >0.70 kU/l) are additionally adjusted for floor 
type of dust sampling. 



with each other (around ;■ = 0.90) and also with amount of 
dust (around ;• = 0.85), but correlations between total micro- 
bial quantity scores and qPCR results were lower (around 
r = 0.30), except for Mycobacterium and the PenAsp group 



(around /• = 0.50). Diversity score had low correlations with 
amount of dust, all microbial markers, and total microbial 
quantity scores (around ;• = 0.20), but higher with Wallemia 
sebi (r = 0.66). 
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In general, the associations of single microbial markers with 
asthma ever and respiratory symptoms were mostly nonsignif- 
icant, while the shape of the association varied (Fig. 1). Total 
microbial quantity score had a significant inverted-U-shaped 
association with asthma ever: the highest risk for asthma ever 
was found with the medium level of total microbial quantity 
score and the lowest risk with the highest quantity score 
(Fig. 2, Table 3). The association was independent of the 
diversity score. The associations were quite similar to current 
asthma (Table 3) and when using other microbial markers 
than ergosterol, endotoxin, and muramic acid to create the 
total microbial quantity score (Table S2) as well as when 
using EPS as a marker of fungi (data not shown). Total 
microbial quantity scores tended to be inversely associated 
with respiratory symptoms (Table S2). The inverted-U-shaped 
association between total microbial quantity score and 
asthma ever tended to be stronger after adjusting for the 
amount of dust (data not shown). Adjustment for dust also 
strengthened the associations between asthma, and Strepto- 
myces, endotoxin, PenAsp group, fi-D-glucan and ergosterol, 
but not the other single microbial markers (data not shown). 

The associations between total microbial quantity score and 
any wheezing were similar among children who were sensitized 
against at least one of the tested inhalant or food allergens at 
the age of 6 years and who were not sensitized (P- value = 0.66 
for interaction term). In contrast, in the stratified analyses, the 
inverted-U-shaped associations between the quantity scores 
and asthma ever were always stronger among nonatopic than 
among atopic children, but the P-values for interaction ranged 
between 0.08 and 0.28 (data not shown). There were no inter- 
actions with birth cohorts (P-values > 0.27 for interaction 
terms). Associations between total microbial quantity score 
and any wheezing were also similar in farms and in rural/sub- 
urban areas (P-value = 0.58 for interaction term). Due to the 

6 - 

5 - 
4 - 



2 - 



1 - 



1 st quinlile 2nd quintile 3rd quintile 4th quintile 5th quintile 
Quantity score 

Figure 2 Risk of asthma in relation to total microbial quantity score 
in quintiles (sum of ergosterol, endotoxin, and muramic acid) (aOR). 
Model is adjusted for study cohort, farming status, gender, mater- 
nal history of allergic diseases (hay fever, atopic dermatitis, and/or 
asthma), smoking during pregnancy, the number of older siblings, 
living area of the home, and diversity score. 
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low numbers of children with asthma ever among farmers, 
interaction tests between total microbial quantity score and 
farming on asthma ever could not be performed. 

Diversity score was associated with a decreased risk of any 
wheezing, and a similar tendency was seen with asthma ever 
and current asthma (Table 3). The association with asthma 
ever was, however, fairly sensitive to different adjustments, 
especially adjustment for farming. The association with any 
wheezing was similar among children from farms and from 
rural/suburban areas as well as children who were sensitized 
against at least one of the tested inhalant or food allergens at 
the age of 6 years and who were not sensitized (P-values 
> 0.4 for interaction terms) (data not shown). 

Most of the associations between single microbial markers 
and atopic dermatitis and sensitization to inhalant allergens 
at the age of 1 (Table S3) and 6 years (Fig. 1) were nonsig- 
nificant. Total microbial quantity score was not associated 
with atopic dermatitis or sensitization to inhalant allergen at 
the age of either 1 year (data not shown) or 6 years 
(Table 3). Diversity score was significantly (inverted-U shape) 
associated with sensitization to inhalant allergen at the age of 
6 years (Table 3), while similar tendency, though nonsignifi- 
cant, was seen with sensitization to inhalant allergen at 
1 year (Table S4). No association was found between diver- 
sity score and atopic dermatitis (Table S4). 

Discussion 

The present study suggests that an index of total quantity of 
environmental microbial exposure may predict asthma better 
than single microbial markers, independently of microbial 
diversity score and amount of dust. The diversity score 
decreased the risk of wheezing and was significantly 
(inverted-U shape) associated with sensitization to inhalant 
allergens. Only a few associations were observed between sin- 
gle microbial markers and the outcomes studied. 

This is the first study using an index of total microbial 
quantity score when studying the risk of asthma. Previously, 
one cross-sectional study found a tendency of inverse associa- 
tion between microbial quantity score (LPSio : o-i6:u, ergos- 
terol, and muramic acid) and allergic sensitization (9). 
Typically, epidemiological studies have used only one micro- 
bial marker as a surrogate for total microbial exposure. 
However, it is clear that such assessments provide an incom- 
plete picture of the total quantity and quality of human 
microbial exposure. Due to high correlations, a major diffi- 
culty of earlier studies has been to separate the effects of sin- 
gle microbial markers from each other and from the effect of 
the amount of dust (5, 20). In the present study, total micro- 
bial quantity score predicted the risk of asthma indepen- 
dently of microbial diversity and amount of dust. 

The reason behind the observed association between total 
microbial quantity score and asthma is unclear. A potential 
explanation is that different components of microbial groups 
or genera share common antigen or molecule structures on 
their cell wall surface, which are detected by pattern recogni- 
tion receptors in innate immunity cells and have a profound 
influences on the immune system (21). If total exposure to 
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such structures is the main causal factor explaining the asso- 
ciation between microbial exposure and lower risk of allergic 
diseases, then total microbial quantity score is likely to better 
reflect this exposure than single microbial markers. Varying 
correlations of single microbial markers with total microbial 
exposure would then explain the often weak and conflicting 
findings from earlier studies. In addition, some studies have 
shown synergistic interactions between different exposing 
agents (22), and an index of total microbial quantity score 
may capture such synergistic effects. Data from deep 
sequencing of microbial exposure (23) will give opportunities 
to explore further the above potential mechanisms. 

Total microbial quantity score had an inverted-U-shaped 
association with asthma. Previous studies have suggested that 
the immune system may react first with enhanced response 
and subsequently with tolerance (24). Associations of similar 
shape have also been reported earlier between allergen levels 
and fj-D-glucan, and risk of asthma, wheezing, or atopic sen- 
sitization (25-27). The shape of the association in the present 
study may also relate to some unknown characteristics of the 
present study, given that some single microbial markers also 
showed an association of a similar shape. Further studies in 
larger cohorts and in different settings are needed to better 
define the shape of the association between total microbial 
quantity score and risk of asthma. 

We found a similar inverse association between measured 
microbial diversity score and wheezing, as has been reported 
with asthma in two recent cross-sectional studies (8). The asso- 
ciation with asthma was nonsignificant, though in the same 
direction. There were differences in the microbiological meth- 
ods used between the present study and previous studies. We 
used qPCR, a method that quantitatively describes the pres- 
ence of specific viable and nonviable organisms of prior chosen 
microbes (28). The previous studies either cultivated viable 
microbes (GABRIELA study) or used single-strand conforma- 
tion polymorphism analyses, which is a qualitative rather than 
quantitative method (PARSIFAL study) (8). We calculated a 
diversity score based on 6 fungal and 2 bacterial qPCR assays 
representing either microbial species, genera, or groups, which 
hence only provides a suggestive and rough estimate of the 
overall microbial diversity. Nonetheless, the results are in line 
with previous studies (8) targeting a broader spectrum of 
microbes. Taken together, the studies suggest that microbial 
diversity score may predict asthma and may be a useful marker 
to describe the quality of microbial exposure. 

The observed associations between studied health out- 
comes and single microbial markers should be interpreted 
with caution, as no adjustments for multiple testing were 
done. The most interesting findings were maybe the protec- 
tive associations in the case of Gram-positive bacteria (Myco- 
bacterium spp., Streptomyces spp., and muramic acid), which 
have been reported earlier for muramic acid (6, 29). 

The main strengths of the present study are the prospective 
birth cohort design with high participation rate (varied 
between 80% and 95% in each follow-up) and a good variety 
of measured markers of exposure to environmental microbes. 
The small size of our study cohort is the main limitation of 
the study, and the results need to be confirmed in further stud- 
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ies. The microbial markers and combinations of those used to 
build total microbial quantity scores in the present study 
do not measure total microbial exposure without error, for 
example using muramic acid together with endotoxin overesti- 
mates the amount of Gram-negative bacteria, because muramic 
acid is found also in Gram-negative bacteria, although to a 
lesser extent (30). These weaknesses do not, however, appear to 
be a major problem, as the different quantity scores used in the 
present study showed very similar associations. 

In conclusion, our indicator for total quantity of microbial 
exposure predicted asthma better than single microbial mark- 
ers. Detected associations were independent of the measured 
microbial diversity and amount of dust. Better indicators of 
total quantity and diversity of microbial exposure are needed 
in studies on the development of asthma. 
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